Scalar Mesons in Lattice QCD 
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We explore whether "fo(600) or cr" exists as a pole in QCD, by a full lattice QCD simulation on 
the 8"^ X 16 lattice using the plaquette action and Wilson fermions. It is shown that a low-mass a 
pole exists owing to the disconnected diagram of the meson propagator. A discussion is given on 
the physical content of the cr. 
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Recently there has been renewed and growing interest 
in scalar mesons 0,0 • One of the most noteworthy devel- 
opments in this field is the re-identification of the (t(600) 
after some twenty years. The results of re-analyses of the 
TT-TT scattering phase shift of the scattering (5) matrix 
strongly suggest the existence of the pole of the a meson 
with / = and J^'-' = 0"*""*" in the s channel as well as 
the p pole in the t channel. In confirming the a pole, 
it was essential to respect chiral symmetry, analyticity, 
unitarity and crossing symmetry in constructing the S 
matrix @, 0, 0| . The significant contributions of the a 
pole were also identified in the decay processes of heavy 
particles such as D— + TTTTTr and T(3S') Ttttf [lllaHH- 
The meson called "/o(600) or a" appeared in the 2002 
edition of PDG 0. It should be noted that a recent 
analysis ^3 clearly shows phase motion consistent with 
resonance behavior. 

As stressed above, the significance of the a pole is in- 
timately related to chiral symmetry in QCD; the a may 
be identified as the chiral partner of the tf in the lin- 
ear representation of chiral SU(2)(g)SU(2) symmetry 
Here, it is noteworthy that the k meson with 7 = 1/2 
is reported to exist with a mass of about 800 MeV 
11 111 El; the K is supposed to constitute the nonet 
scalar states of chiral SU(3)(8)SU(3) symmetry together 
with the a. See 'Meson Particle Listings /o(600)' in Ref. 
(pp.01000 1-450 - 453), and references therein. 

Even with the confirmed existence of the low lying 
scalar mesons, there are long standing controversies on 
their nature Are they usual qq mesons, four quark 
states like qqqq ^^li i^^) molecules, rather than 

pre-existing resonances? How large is the mixing among 
them or with glueballs ^B]"^ Although not mentioned in 
0, the a may be a collective qq state described as a su- 
perposition of many atomic qq states [l6|,|,17]: In fact, the 
dynamical models of chiral symmetry breaking, such as 
the Nambu-Jona-Lasinio model 17] and the Schwinger- 
Dyson approaches [ll| , describe the a as well as the tt as 
collective states. 



Confronting these problems with the cr, it would be 
intriguing to explore whether QCD accommodates the a 
and other possible low- lying scalar mesons, as well as to 
examine their nature by a first-principles calculation of 
QCD. 

DeTar and Kogut were the first who measured the 
screening mass of the cr meson in lattice QCD, although 
in the quenched approximation. In the same approxi- 
mation, Alford and Jaffe explored possible light scalar 
mesons described as q^if states 20], while the masses 
of the qq states and their mixing with glueball states 
in the scalar channel were investigated by Lee and Wein- 
garten ^3] • It was pointed out [23] that the so-called dis- 
connected diagram, which is not properly treated in the 
quenched approximation, gives almost the same amount 
of contribution to the propagator as the connected dia- 
gram does; this indicates that using the dynamical quarks 
is essential to study the cr on a lattice. 

A simulation with dynamical quarks was carried out 
by McNeile and Michael [2I], who computed the masses 
of the iso-singlet scalar states given by a superposition of 
qq and glueball states. In their simulation, the cr meson 
was found to be lighter than the tt; hence the relevance 
of the simulation in reality is obscure. 

Since the existence and the nature of the a may be 
related to chiral symmetry and its dynamical breaking, 
it is desirable to examine the behavior of the chiral limit 
of the calculated quantities for a proper description of 
the cr meson. Of course, it is preferable to adopt lat- 
tice fermions which have better chiral properties. Us- 
ing the domain wall fermions, RBC (Rikcn-Brookhavcn- 
Columbia) collaboration examined the cr propagator on 
the lattice, unfortunately, however, in the quenched ap- 
proximation and with an approximate estimate of the 
disconnected diagram based on the chiral perturbation 
theory 

In this letter, we present a lattice calculation of the 
scalar particles using full QCD with the inclusion of the 
disconnected diagrams. We employ the most standard 
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lattice QCD numerical techniques, and show that there 
exists a o pole, whose mass is similar to that of the p 
meson for the quark masses used for the simulation but 
tends to become as low as IttIt^ in the chiral limit. We 
shall show that the incorporation of the disconnected di- 
agram is responsible for realizing such a low mass of the 
a meson. We shall briefly discuss the physical content of 
the a meson. We shall also mention that our simulation 
suggests the existence of the k meson, another member 
of the nonet of the scalar mesons. 

We adopt the following interpolation operator for cre- 
ating the (T meson with / = and J^*^ = 

-i-) = Y.Y. 71 ' w 

c— 1 a— 1 

where u (d) denotes the up-quark (down-quark) operator 
with c and a being the color and Dirac-spinor indices, 
respectively. Then the a meson propagator reads 

Giy,x) = {Ta{y)a{xy) 

= ^ J DUDuDuDdDd 

X a{y)a{x)'^ ^-So-uWu-dWd ^ (2) 

Here W~^^s are the u- and d-quark propagators, C/'s the 
link variables of the gluon, and Sg the gauge action. By 
integration over the quark fields, the a meson propagator 
is reduced to 

G(y,x) - -{TvW-\x,y)W-\y,x)) 

+2{{a{y)~{a{yma{x)-{a{x))))^ (3) 

where 

a{x) = TT:W-^{x,x), (4) 

with "Tr" being the trace operation over the color and 
Dirac indices. Here we set the u and d quark propagators 
to be the same. Equation ||2Jl shows that the a propagator 
consists of two terms: one corresponds to a connected 
diagram and the other to a disconnected diagram. Since 
the vacuum expectation value {(j{x)) does not vanish, it 
should be subtracted from the a operator. We employ 
Wilson fermions and the plaquette gauge action. We 
perform a full QCD simulation in which the disconnected 
diagrams are included. 

As for the simulation parameters, we first note that 
CP-PACS performed a fuU-QCD calculation of the hght 
meson spectroscopy with great success Therefore, 
we use the same values of the simulation parameters as 
those used by CP-PACS, i.e., /3 = 4.8 and the hopping 
parameter, k = 0.1846, 0.1874 and 0.1891, except for the 
lattice size; our lattice size is 8'^ x 16. We shall use the 
point source and sink, which together with the smaller 
lattice size, leads to larger masses due to a mixture of 



higher mass states. In other words, the masses to be 
obtained in our simulation should be considered as up- 
per limits. We show in Table 1 the value of m^/mp 
together with that of CP-PACS. They are consistent 
within the error bars. We generate the gauge con- 



Table 1: Summary of the results. 



K 


0.1846 


0.1874 


0.1891 


statistics 


1110 


860 


730 


niTy/nip 


0.8291(12) 


0.7715(17) 


0.7026(32) 


m,r/mp ^> 


0.825(2) 


0.757(2) 


0.693(3) 




1.6(1) 


1.34(8) 


1.11(6) 




2.40(2) 


2.44(3) 


2.48(4) 



'Number of configurations separated by 10 trajectories to 
each other. ^'CP-PACS. ^'Our result, mconncct is the a mass 
estimated only from the connected part. 

figurations in full QCD using the hybrid Monte Carlo 
(HMC) algorithm. The first 1500 trajectories are up- 
dated in quenched QCD, then we switch to a simulation 
with the dynamical fermion. The next 2000 trajectories 
of HMC are discarded for thermalization and the ct, tt 
and p propagators are calculated every ten trajectories. 

It is not an easy task to evaluate the disconnected part 
of the propagator, since one must calculate TrW^~^(x,x) 
for all lattice sites x. We used the Z2 noise method to 
calculate the disconnected diagrams, <j(x), and the sub- 
traction terms of the vacuum (cr) . Each of these terms is 
the order of ten, and ((a — {(t)){(t — {cr))) becomes less 
than 10~^, as shown in Fig^ Therefore, in order to ob- 
tain the signals correctly as the difference between these 
terms, numerical simulations with a high precision and 
careful analyses are necessary. One thousand random Z2 
numbers are generated for each configuration: We refer to 
our previous report |2^ for the relationship between the 
amount of Z2 noise and the achieved accuracy. Gauge 
configurations are created by HMC in a vector super- 
computer, while most of the disconnected propagator is 
calculated in a parallel machine. 

The propagators of the tt, p and cr for k = 0.1874 
are shown in Fig[21 From our results, we estimate the 
critical value of the hopping parameter and the lattice 
space as Kc = 0.195(3) and a = 0.207(9) fm, respec- 
tively, which are to be compared with the CP-PACS val- 
ues, Kc = 0.19286(14) and a = 0.197(2) fm. The mass 
ratios of rriTr/mp and ma- /nip together with mconnect/wp 
are summarized in Table 1, where TOconncct is the mass 
evaluated only from the connected diagram. We have 
extracted the tt and p masses using two-pole formula 
from their propagators. The small errors of the tt and 
p propagators indicate the high precision of our simula- 
tion. Our result for the mass ratios 771^^/ nip in Table 1 
is in a good agreement with that of CP-PACS, with un- 
avoidable small differences owing to the different lattice 
sizes. To extract the a meson mass, we have used the 
one-pole ansatz, since its propagator has large error bars 
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at large t ^ 6 in spite of our high statistics; it imphes 
that the obtained value gives an upper limit of the scalar 
meson mass. Figure Eland Table 1 show that the mass of 
the a is of the same order as the mass of the p; we remark 
that the 27r threshold is higher than the a in our simu- 
lation. We have checked that the effective mass method 
gives essentially the same result as that given above. 
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FIG. 1: Propagators of the connected and disconnected dia- 
grams of the a for K = 0.1874. 



with the fact that the oq is heavier than the p. Figure Q] 
shows that the disconnected part dominates the a prop- 
agator. By comparing Fig|21with Fig^ we see that the 
(T as a light meson results from the disconnected part 
of the a propagator with the background vacuum con- 
densate subtracted, which are odd with the constituent 
quark model. In the naive sense, the connected quark di- 
agram corresponds to the constituent quarks in the SU(3) 
non-relativistic quark model where the OZI rule is satis- 
fied. This may give a clue to clarify the physical content 
of the a meson, as will be discussed later. The mass of 
the connected a (the qq) shown in Table 1 is almost 2.5 
times of the p mass and exhibits only a weak dependence 
on the hopping parameter, suggesting an irrelevance of 
chiral symmetry to the ag meson. 
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FIG. 2: Propagators of the n, p and cr for k = 0.1874. 



Figure ^ shows the individual contributions of the con- 
nected and disconnected parts of the a propagator, which 
reaches a plateau ai t > 5, since the precision of our cal- 
culation is hmited to 0{G{t)) ~ 10"'*. The connected 
part shows clear signals of a rapid damping with small 
error bars. We note that the connected part of the a 
propagator can be regarded as the one of the ao meson 
provided that the difference between the M-quark and the 
d-quark is neglected. Therefore, the rapid damping of 
the connected part of the a propagator is in accordance 




FIG. 3: m^, nip, rria and 2m,r in the lattice unit as a function 
of the inverse hopping parameter. The chiral limit is given by 
Kc = 0.195(3). 



In Fig|2| we display m^, rrip, iria- and 2m^ in lattice 
units as a function of the inverse hopping parameter. As 
the chiral limit is approached, the a meson mass obtained 
from the full a propagator decreases and eventually be- 
comes smaller than the p meson mass in the chiral limit. 

We have also calculated the k meson using the same 
configurations. We take the same hopping parameter val- 
ues for u and d quarks and use the following three values 
for the s quark hopping parameter, Kg = 0.1835, 0.1840 
and 0.1845. The ratios rrii^/mK* and rriK/mK' at the 
chiral limit are listed in Table 2. One may notice that all 
the three values of Kg fairly reproduce the physical mass 
ratio, mK/mK* in the chiral limit of the u- and d-quarks. 
A more notable point is that the resultant k mass hardly 
changes with the variation of Kg and is almost twice as 
heavy as the K*; i.e., the k is not degenerated with the a, 
although they should belong to the same nonet. The ori- 
gin to lift the octet-singlet degeneracy may be attributed 
the following facts; because of the strangeness content. 
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the K propagator is solely composed of a connected di- 
agram and contains no disconnected part, the latter of 
which was the origin of the light mass of the a. In fact, 
the K mass listed in Table 2 is almost the same as the Qq 
mass estimated using the connected part of the a with the 
corresponding hopping parameters. Thus we can under- 
stand why the a is light while the k is heavy, consistently. 

Table 2: The mass ratios mi^/rriK* versus mK/mK' at 
Kc = 0.1945(29). The s quark hopping parameters were 
taken to be 0.1835, 0.1840 and 0.1845. 



(No. 11440080, No. 12554008, No. 12640263 and No. 
14540263) and ERI of Tokuyama Univ and DOE grants 
DE-FG02-96ER40495. Simulations were performed on 
SR8000 at IMC, Hiroshima Univ., SX5 at RCNP, Osaka 
Univ., and SR8000 at KEK. 



s quark hopping parameter 


0.1835 


0.1840 


0.1845 


rriK/mK* 


0.639(6) 


0.631(6) 


0.623(6) 


rriK /rriK- 


2.039(43) 


2.037(43) 


2.044(44) 



The dominance of the disconnected diagram in the a 
meson propagator in contrast to the octet scalar mesons 
suggests possible physical contents of the light a: The 
disconnected diagram includes the process qq ^ G 
where G denotes a glueball states. This might suggest 
an importance of the glueball mixing in the light a if 
the glueball states were not heavy. Instead, we notice 
that such a fluctuation process of the qq states with the 
heavy G giving an effective vertex for qq — s- qq may form 
a collective state mentioned before. One may also notice 
that the disconnected diagram contains four quark states 
composed of diquark and anti-diquark states, i.e., qq-qq. 
Clearly more works are needed to elucidate the physical 
content of the scalar mesons. 

We have reported the first study of the scalar mesons 
based on the full QCD lattice simulation with dynami- 
cal fermions, including analysis of the disconnected di- 
agram effects. We have used the most standard lattice 
QCD techniques, which have worked well for the estab- 
lished mesons and baryons, and clarified the accuracy and 
statistics required to obtain signals in the scalar channel. 

Our results indicate the existence of a light isoscalar 
J^'-^' = 0++ scalar meson, i.e., the a meson with a mass 
of almost the same order as that of the p meson (see Table 
1). We have also shown that the k meson belonging to 
the flavor octet is almost twice as heavy as the K*. 

Of course, the results reported here are in far from 
final quantitative stage because our quark masses are 
much heavier than those in the real world and our scalar 
mesons cannot decay on our lattice. Nevertheless, it 
is now clear that the a meson as well as other scalar 
mesons can be studied in the lattice QCD. We hope that 
the present work prompts other lattice studies which will 
give a deeper and more quantitative understanding of the 
scalar mesons, especially the a. 
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